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QED in 2+1 dimensions has long been studied as a model field theory which exhibits both asymptotic freedom 
and non-trivial IR behaviour. There is also a trend towards viewing it as a candidate low energy efi'ective theory 
for the pseudogap phase of high temperature superconductors. One feature of these theories is their lack of 
isotropy in the x and y directions (a common feature of Dirac theories in condensed matter systems). A model 
motivated by this work is outlined, and tentative evidence presented for chiral symmetry restoration as the relative 
anisotropy is increased. 



1. Introduction 

One of the outstanding questions in the theory 
of condensed matter systems is the nature of the 
so-called pseudogap phase of the cuprate family 
of high temperature superconductors. 

This phase, which is reached from the super- 
conducting (dSC) phase by a combination of in- 
creased temperature and reduced doping (x), is 
characterised by the persistence of an energy gap 
but without the corresponding sharp quasiparti- 
cle peak — that is, an energy gap but no su- 
perconductivity. Traditionally, approaches to this 
problem begin with a model describing the anti- 
ferromagnetic (AF) phase as a starting point and 
attempt to derive the properties of the pseudo- 
gap from it; however, a more recent trend, due to 
P, begins from the dSC phase instead. Inspired 
by this approach, the authors of |2] and |2] have 
derived a low temperature effective theory which 
very much resembles quantum electrodynamics in 
2+1 dimensions (QED3) and this in turn has mo- 
tivated our research into the area. 

2. From dSC to QED3: A brief review 

Following J2, we begin with the action of a d- 
wave superconductor at T 0: 

S^T ^ [(iLj„ - ^k)4(k,w„)cCT(k,a;„) 
- |A(k)ct(k,c^„)c^,(-k,-^„) 



+ h.c. + 0{c^)] (1) 

where c\,, Ca are the operators, cr = +, — signi- 
fies spin, ujn is the fermionic Matsubara frequency 
and O(c^) are terms governing short range inter- 
actions. A(k) is the gap function, which has d- 
wave symmetry. 

We expand around the four nodes of the 
d-wave Fermi surface, and rewrite in terms of the 
fields *|(q,(.j„) = (c^k,a;„),c_(-k, -Lj„),c^k - 
Qi,tj„), c_^(-k + Qi, -tj„)). Here, i = 1,2 and 
the wave vectors linking the nodes into diagonal 
pairs are Qi = 2Kj, where is the distance of 
the node from the origin in momentum space. 

For k = Ki + q, |q| <C |Ki| we can linearise cer- 
tain parameters as = vpq-j^ + 0{q^), A(k) = 
VAQy + 0{q^). vp, VA are referred to as the 
Fermi and Gap velocities respectively. If we de- 
fine S — ^yvpvA and the anisotropy parameter^ 
K = ^ , as is done in 0j , we can make the replace- 
ments vfPx ^ "T^Pa: ^^"i VAPy &\fK.Px vo. the lin- 
earisation; this makes the comparison of results 
with the predictions within 0] easier. Transform- 
ing the result of our manipulations into position 
space we have: 



+ (1^2,x^y)} 



(2) 



^This parameter depends on the doping x of the super- 
conductor. HI! states that the value of this parameter at 
optimum doping is about 10. 



1 



2 



This describes a noninteracting system; the phase 
degree of freedom in the order parameter A(k) al- 
lows us to add gauge fields and interaction terms. 
Detailed workings out of this are given in \'2l'6\ . 
but for our purposes, the more heurisitic accounts 
given in |5I6| will suffice. They proceed as follows: 
the phase d.o.f. allows the formation of vortices. 
These vortices, if they are allowed to condense, 
will disorder the dSC phase, placing us within the 
pseudogap. The photon in 2+1 dimensions has 
the same physical d.o.f. as the Goldstone boson 
produced by vortex condensation, so it is natu- 
ral to model the effects of the phase vortices as 
a minimally coupled U(l) gauge theory^. Doing 
this, we can write the anisotropic QED3 action 
of |2I3I4| ■ where a^j. is the photon field, and e is 
expectation value of the vortex condensate ($): 



S = / d'^rdT{^i[jo{dr+tar) 



--li{dx 



Sy/Kj2idy + iay)]^i + {1 ^ 2,x ^ y) 
1 



(3) 



3. Relevance of chiral symmetry breaking 

One of the more interesting questions regarding 
the phase diagram at T = is whether the pseu- 
dogap exists as an intermediate phase between 
the AF phase and the dSC phase. j3] and T show 
that the chiral condensate (^^E") can be identified 
with a Spin Density Wave order parameter (which 
becomes a full AF phase as x decreases). Conse- 
quently, = signals the existence of the 
pseudogap. Note that at no point can our model 
enter the dSC phase; the addition of the gauge 
fields to our action presupposed that we had left 
it. 

The much vexed question of the critical num- 
ber of flavours (Nc) required for chiral symme- 
try breaking in QED3 now becomes relevant. If 
Nc < 2 Vk, then no pseudogap exists as an in- 
termediate phase at T = 0. If A'^c > 2 Vk, then 
it does. Schwinger-Dyson techniques have given 



^The correct approach for modelling this on the lattice is 
to use noncompact formulations of U{1). Due to instan- 
ton formation, compact formulations do not preserve flux 
symmetry ^ 0), hence the photon in such a model 

is massive 151. 151. 
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Figure 1. (^'5') versus k for mo between 0.01 and 
0.05. 

values for A^c of 00, || and ~ 4 (see references 
cited in Uj). A thermodynamic argument due to 
[H] has A^c < |, whilst |7j has an argument that 
A'c > 4 due to a correspondence with the Thirring 
model. Lattice calculations [7] have failed to shed 
much light on the issue; chiral condensates (after 
the appropriate limits are taken, and allowing for 
the finite size effects that bedevil QED3 simula- 
tions) are very small, which indicates either sym- 
metry restoration, or symmetry breaking with a 
very small condensate. It is hoped that despite 
these issues, simulations of the anisotropic case 
might allow general conclusions to be drawn as 
to the nature of the phase diagram. 

4. The Lattice Model 

We simulate the following staggered fermion 
action using a Hybrid Monte-Carlo algorithm on 
a 16"^ lattice: 
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The fermion matrix is defined as follows: 

3 

Mx,x' = E ^M(2^)[^a;',x4-A^2;M " ^.c] + TOqI^^i. (5) 

Here, ^^(a;) = X^t]^{x). The anisotropics enter 
as Xfj,: Ai = k"2, A2 = k^, A3 = 1. Vt^i^) = 
(^_l^xi+...+x^, is the Kawomoto-Smit phase. Uxii = 
exp{i6x^i), and the noncompact action is defined 
by e^,(x) = A+9,ix) - A+9^ix). N = 1, which 
gives us Nf = 2 in the continuum limit as re- 
quired. Note that this lattice action is intended 



3 



mo = 0.04 --X-J :■ 
1 - nio = 0.03 r -y^ 

mo = 0.02 :- □ -i ■ ' 

Q g _ mo = 0.01 - . ^] 

|0.6 - * + »^ ■ 

0.4sr i i i ^ H S ^ 

[] H □ E3 H ^ 

0.2"- 

g I I I I I I I I I 

123456789 10 

K 

Figure 2. m.^. versus k for nio between 0.01 and 
0.05. 

to correspond to the anisotropy structure of the 
first flavour in with 6=1. 

5. Preliminary results and comments 

Here we present preliminary resuhs concern- 
ing the chiral condensate, pion mass and renor- 
malised anisotropy k^- 

The chiral condensate, shown in Figure ^ 
exhibits a marked decrease as k, the bare 
anisotropy, increases. This could be either a 
phase transition at k « 5 or a crossover; work 
needs to be done regarding the susceptibihties in 
order to clarify this. 

The pion mass, extracted from the timeslice 
propagator, is shown in Figure |21 We must cau- 
tion that the error estimates are somewhat con- 
jectural; the pions appear to be too light for 
the lattice. Naively, one would expect that the 
pion mass would be unaffected by the change in 
anisotropy, as k does not enter lO in the r di- 
rection. However, as we see, the mass increases 
roughly quadratically with respect to k. This has 
yet to be explained; though it is plausibly due to 
some finite size effect, genuine physics can't be 
ruled out at this stage. 

Lastly, we present the renormalised anisotropy 
Kr (Figure OJ), which is defined by the ratio of 
the effective pion masses (the inverted correlation 
lengths) in the x and y directions. j3] present re- 
sults from a large-A'^ expansion suggesting that 
K is an irrelevant parameter. However, here we 
can see that Kr is greater than k for 1 < k < 3, 
which suggests that it is a relevant parameter. 
Whether this contradicts 4 or not is not cer- 
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Figure 3. Kr versus k for mo between 0.04 and 

0. 05. 

tain — q.v. the discussion at the end of the pre- 
vious section on the anisotropy structure of our 
model. As has been emphasised, the results given 
here are very preliminary; certain immediate re- 
search directions have been outlined in the dis- 
cussion above. It will be necessary to simulate on 
larger lattices with different values of (3 (how Kr 
depends on this parameter is an intriguing ques- 
tion) in order to take the continuum limit, a task 
that has often proven difficult with QED3. We 
hope, however, that we can still shed some light 
on the behaviour of this novel field theory and its 
connection with a fascinating area of condensed 
matter physics. 
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